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EXTRACTION FROM FLIGHT DATA OF LONGITUDINAL AERODYNAMIC 
COEFFICIENTS FOR F-8 AIRCRAFT WITH SUPERCRITICAL WING 

By James L. Williams and William T. Suit 
Langley Research Center 

SUMMARY 

A parameter -extraction algorithm was used to determine the longitudinal aerody- 
namic derivatives from flight data for the F-8 aircraft with supercritical wing. The flight 
data were the responses to horizontal-tail pulses and the data used were for Mach num- 
bers of 0.80, 0.90, and 0.98. 

Results of this study showed that a set of derivatives were determined which yielded 
a calculated aircraft response almost identical with the measured response in flight. 

In addition, results of this study showed that the trends of the converged derivatives 
with Mach numbers were generally similar to the trends obtained by using a Newton- 
Raphson method. At the highest Mach number, the converged damping derivative was 
substantially lower than the value of the damping derivative obtained from preliminary 
tests in the Langley 8-foot transonic pressure tunnel. This lack of agreement probably 
resulted from the difference in Reynolds number between the flight and wind-tunnel tests. 

INTRODUCTION 

Analytical and simulator studies of the flight and handling qualities of aircraft 
require that accurate estimates of the aerodynamic parameters be used if the results are 
to be valid. To provide aerodynamic parameters for analytical and simulator studies and 
to also provide numerical values for comparison with wind-tunnel data and theoretical 
estimates, parameters have been extracted from flight data for many yeais. Results 
from recent studies made at the Langley Research Center relating to parameter extrac- 
tion for several aircraft are reported in references 1 to 4. 

Various techniques for parameter extraction are presented in reference 5. The 
technique and program used in extracting the parameters in this study are those of the 
maximum likelihood method of reference 6. 

The purpose of the present study is to extract longitudinal aerodynamic parameters 
from flight tests of the F-8 supercritical-wing aircraft at Mach numbers of 0.80, 0.90, 
and 0.98. 
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SYMBOLS 


Ax>Az linear accelerations, g units 

Cm pitching-moment coefficient about Y body axis 

r - 9C m 
v ' m q “ 

c m fl = P er radian 





8C[j] 

— ~ > per radian 

95 h 


C X force coefficient along X body axis 

OCx 

per radian 


C Z force coefficient, along Z body axis 

_ 9C Z 

c Z a = per radian 



0C Z .. 

= ■' P er radian 

06 h 


c mean geometric chord, m 

DE = 6 h - 5h,t» rad 

g gravitational acceleration, lg = 9.81 m/sec 2 

h p pressure altitude, m 
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I X» I Y» I Z 

J xz 

*h 

N 

P 

q 

**=k 

R _ 1 

U ~W^ 

Rw = — 
^ W w 

r 

S 

u 

V 

v 

W 

Wq,W u ,W w 

w 

X 


moment of inertia about the roll, pitch, and yaw axis, kg-m 2 
product of inertia, kg-m 2 

distance from aircraft, center of gravity to quarter-chord point of mean 
aerodynamic chord of horizontal tail, m 

number of data points 

roll rate, rad/sec 

pitch rate, rad/sec 


yaw rate, rad/sec 
wing area, m 2 
velocity component along 
resultant velocity, m/sec 


velocity component along 
vector of aircraft states 


X body axis, m/sec 


body axis, m/sec 


w state variables 
body axis, m/sec 


velocity component along Y 
aircraft weight, N 
weighting value of q, u, and 

Z 
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a 


angle of attack, deg or rad 


0 angle of sideslip, deg or rad 

6jj horizontal-tail deflection (positive for trailing edge down), deg or rad 

6 pitch attitude angle, deg or rad 

p mass density of air, kg/rn^ 

<p roll attitude, rad 

Subscripts: 

c computed 

h horizontal tail 

m measured in flight 

t indicates state at trim conditions 

The dot over a symbol denotes rate of change with respect to time. 

EQUATIONS OF MOTION 

The equations of motion used in this study are referred to a body-axis system 
(fig. 1) and are as follows: 

X Translational direction: 

i = rv - qw - g sin S ♦ I p y2s | fc x>t ♦ C Xq , (a - a t ) 

L_ _ 

Z Translational direction: 

w = qu - pv + g cos « cos <f> + Ipv2 S |fc 2jt + C Za [a - « t ) + Cz 5h (6 h - 5 hjt ) 
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Pitching moment: 


q = pr + ^p(r 2 - p2) + |^l“[^m 3 t + C m a (“ ' a t) + C m q §| 

+ Cm a + Cm 6h( 5h " 6h,t )J 

The lateral variables v, p, and r were found to be small, and hence the nonlinear 
terms in these equations were negligible for these calculations. 

DESCRIPTION OF AIRPLANE 

The airplane for which the aerodynamic parameters were extracted is the F-8 air- 
plane on which the original wing was replaced with a supercritical wing. The aircraft is 
a single-seat high performance airplane with a single jet engine embedded in the fuselage 
and a unit horizontal tail. Figure 2 is a photograph of the airplane. Pertinent geometric 
characteristics of the airplane are given in table I. 

Flight instrumentation appropriate to this study included the following items: 

(1) Pitch-rate gyro 

(2) Angle -of -attack indicator 

(3) Altimeter 

(4) Total velocity indicator 

(5) Horizontal-tail position indicator 

(6) Magnetic tape recorders 

(7) Accelerometers 

The full-scale range of the flight instruments and their accuracies are given in table II. 

FLIGHT TESTS 

The data which are used in the present investigation were obtained from in-flight 
measurement of the airplane response to a horizontal-tail input. These flights were 
made at the NASA Flight Research Center as part of a general evaluation program of the 
F-8 supercritical-wing aerodynamics. Data were obtained at Mach numbers of about 0.80, 
0.90, and 0.98. The control input used to generate the longitudinal motion were one or 
more horizontal -tail pulses. All data used in the study were reduced at Flight Research 
Center and have been corrected for bias and displacement of the measuring instrument 


5 



from the aircraft center of gravity. A list of the test conditions and mass characteristics 
are presented in table III. 

DERIVATIVE -EXTRACTION PROCEDURE 

The parameter-extraction procedure used in this study is an iterative technique 
which utilizes the maximum likelihood method to estimate the stability and control 
parameters. This method uses the likelihood function which, when maximized, provides 
the following information: 

(1) The parameter changes which are used to update the parameter 

(2) The covariance matrices whose elements are proportional to the estimated 

standard deviations and the pairwise correlation coefficients for the param- 
eters and the states 

(3) The performance index function J which is an indicator of the fit between 

measured and calculated motions 

Details of the method are given in reference 6. 

The iterative technique produces a set of estimated derivatives which, when used in 
the equations of motion, provide the best fit to the time variation of the aircraft motion 
measured in flight. The criterion for the best fit is the performance index J which is 
defined as: 


N 


J = det 


iZ( 

L i=l 


x i,m " x i,c) ( X i,m “ x i,c) 


where det means determinant, T means transpose, and X is the vector describing 
the state of the aircraft. Generally, the performance index J becomes smaller with 
successive iterations. The iteration procedure is stopped when the value of J does not 
change appreciably for several successive iterations. The components of the vector X 
are the state variables u, w, q, 0 , and A£. The linear acceleration Ax was not 
included in the X vector since this variable was generally small. The quantities u m 
and w m were not measured directly but were obtained from the measured total velocity 
and angle of attack through the use of equations: 


u m = y m cos “m 

w m = v m sin or m 
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Maximization of the likelihood function yields the covariance matrix for the meas- 
urement noise based on the current nominal solution- (See ref. 6.) This matrix gives the 
variances (or standard deviations) of the differences of the measured state and the nominal 
solution. The inverse of this matrix is the weighting matrix used in the parameter- change 
equations. In this investigation, the diagonal form of the weighting matrix was used and 
the diagonal elements can be expressed as the squares of the difference between the meas- 
ured and calculated data. For example, the weight for the state variable u is expressed 
as 


X ( Um ~ Uc )i 2 

"u i=l 

Similar equations are obtained for l/R^ and l/Rq^. 

Initial values of the state variable were obtained from the flight records for the time 
period just prior to a control input. Initial values of the aerodynamic derivatives were 
obtained from a preliminary investigation in the Langley 8 -foot transonic pressure tunnel 
and these values are listed in table IV. Although the tabulated derivatives were for a Mach 
number of 0.8, they were used as initial values for all the work presented herein. 

RESULTS AND DISCUSSION 

The aerodynamic derivatives of table IV were used as starting values in the 
parameter-estimation process for each flight condition. Preliminary results indicated a 
high correlation between Cz g and and also between C mq and C m ^. These 

correlations had been observed in several other parameter-extraction studies and, there- 
fore, were expected. The correlation problem was circumvented by fixing C m ^ at -7.00 
(the estimated value from table IV) and letting the computer search for C niq . The corre- 
lation between C mgh and Cz 5 h was circumvented by using the geometric relationship 

Cz «h = 4 Cm6 h 

and searching for C ms . 

°h 

The preliminary results also showed that the derivative Cx a had a very large 
variance, which meant that it was not well defined. It was also found that Cx a could be 
varied over a rather large range of values without appreciably affecting the computed 
time histories or the other extracted derivatives. The inaccuracy in this derivative 
resulted from the fact that only small changes in the forward velocity u of the aircraft 
occurred after the control input (Ax small) and this derivative occurs in the equation 
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for forward velocity. The strategy employed in selecting a value for Cx a was to try a 
range of values, and to select that value which resulted in the smallest variances in the 
states at convergence. The value selected was Cx a = 0.58. It appeared to be of reason- 
able magnitude, and was therefore held constant during the iterative process for each 
flight condition. 

Calculated time histories obtained after convergence, and with the constraints noted, 
are shown in figures 3, 4, and 5 for Mach numbers of 0.80, 0.90, and 0.98, respectively. 
The calculated time histories are almost identical with the flight records. 

The extracted derivatives and their estimated standard deviations are listed in 
table V. The extracted derivatives appear to be well defined. There are no standard 
deviations listed for Cx ft , Cz Q , and C m ^ since these were constrained to the tabu- 
lated values, as mentioned previously. 

The more important extracted derivatives of table V are presented in figure 6 as 
functions of Mach number. Also shown are values from preliminary tests in the Langley 
8 -foot transonic pressure tunnel and extracted derivatives from flight tests calculated by 
using a Newton- Raphson method (ref. 7). In general, the extracted derivatives (the 
present study and the Newton- Raphson method (ref. 7}) are in good agreement with the 
wind-tunnel values. However, the extracted effective damping-in-pitch parameter 
" ( c mq + c m^) has a trend with Mach number which is opposite (lower value) that 
obtained from wind-tunnel tests in the Mach number range from 0.90 to 1.00. One reason 
for this difference can be attributed to a Reynolds number effect; since, in general, the 
Reynolds number for the present tests was about 10 times the Reynolds number of the 
wind-tunnel tests at these Mach numbers. 

CONCLUDING REMARKS 

£U=*} 

A parameter-extraction algorithm was used to determine the longitudinal aerody- 
namic derivatives from flight data for the F-8 aircraft with supercritical wing. The 
flight data were the responses to horizontal-tail pulses and the data used were for Mach 
numbers of 0.80, 0.90, and 0.98. 

Results of this study showed that a set of derivatives were determined which yielded 
a calculated aircraft response almost identical with the measured response in flight. 

In addition, results of this study showed that the trends of the converged derivatives 
with Mach number were generally similar to the trends obtained with the Newton-Raphson 
method. At the highest Mach number, the converged damping derivative was substantially 
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lower than the value of the damping derivative obtained from preliminary tests in the 
Langley 8-foot transonic pressure tunnel. This discrepancy appears to be caused by a 
difference in Reynolds numbers between the flight and wind-tunnel tests. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 10, 1974. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF F-8 AIRCRAFT 


Fuselage: 

Length, 16.09 

Wing: 

Area, m2 25.50 

Aspect ratio 6.77 

Span, m 13.14 

Mean geometric chord, m 2.08 

Vertical tail: 

Area, m2 10,13 

Aspect ratio 1.5 

Span, m 3.89 

Rudder: 

Area, m2 1.17 

Horizontal tail: 

Area, m2 8.68 

Aspect ratio 3.5 

Span, m 5.52 

Tail length, center of gravity to quarter -chord point 

of mean geometric chord, m 5.31 
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TABLE II.- INSTRUMENT RANGES AND ACCURACIES 


Instrument 

Range 

Accuracy 

Pitch rate 
Angle of attack 
Pressure altitude 
Total velocity 
Horizontal-tail position 
Normal acceleration at c.g. 
Longitudinal acceleration at c.g. 
j Pitch angle 

±40 deg/sec 
-5° to +25° 

0 to 18 000 m 
0 to 360 m/sec 
-25° to 70 
-lg to +4g 
±0.5g 
±30° 

±0.8 deg/sec 
±0.3° 

±0.024 kN/m 2 
±0.7 m/sec 
±0.4° 

±0.06g 

±0.01g 

±0.6° 


TABLE III.- TEST CONDITIONS AND MASS CHARACTERISTICS 


Weight, 

N 

Mass, 

kg 

ix> 

kg-m^ 

ly, 

kg-m^ 

l Z> 

kg -m 2 


l XZ> 

kg-m 2 

m 

*h> 

m 



Mach number, 0.80 


105 981 

10807.0 

20 519.9 

125574.1 

139 580.1 

4559.5 

11264.0 

5.31 



M 

[ach number, 0.90 


104 722 

10 678.7 

20511.2 

125 305.9 

139319.8 

4514.6 

11274.2 

5.30 


Mach number, 0.98 



101 982 

10 399.3 

20492.3 

124 660.9 

138 693.6 

4406.5 

13 856.0 

5.27 


TABLE IV.- STARTING VALUES OF AERODYNAMIC DERIVATIVES 


CXq, 

c Z a 

Cz «h 

Cm « 


C mq 

° m ®h 

-0.29 

-6.13 

-0.91 

-1.14 

-7.0 

-41.6 

-2.69 









TABLE V.- EXTRACTED AERODYNAMIC PARAMETERS 
AND STANDARD DEVIATION 

Parameters and standard deviations a 


Coefficients 


at Mach number of — 




0.81 

0.90 


0.98 

c X,t • ■ • • 


0.03 

(0.00) 

0.02 

(0.00) 


0.02 

(0.00) 

C Xq! ‘ • 


^0. 58 


b 0.58 



b 0.58 


c z,t • • • • 


-0.44 

(0.00) 

-0.44 

(0.00) 


-0.44 

(0.00) 

c Z a • • • • 


-6.02 

(0.08) 

-6.95 

(0.13) 


-7.49 

(0.24) 

Cz V ■ • ■ 


b -0,92 


b -1.01 


b 

-0.97 


Cm,t • • • ■ 


0.003 

(0.00) 

0.003 

(0.03) 


0.00 

(0.00) 

Cm a • • ■ • 


-1.62 

(0.01) 

-2.04 

(0.10) 


-3.31 

(0.02) 

Cm d? ■ • • • 


b -7.00 


b -7.00 


b 

-7.00 


C mq .... 


-27.82 

(1-04) 

-29.17 

(1.83) 

- 

■16.59 

(2.60) 

c m 6h . • ■ • 


-2.34 

(0.02) 

-2.58 

(0.73) 


-2.46 

(0.07) 


1 — — i 

a Standard deviations are given in parentheses. 
b Used as constants. 



Y 




Figure 1.- System of axes. Positive directions of forces, moments, 
and angles are indicated by arrows. 
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Figure 2.- Photograph of aircraft for which derivatives were extracted. 





















Figure 4.- Comparison of flight data with time histories 
computed by using the aerodynamic parameters of 
table IV for elevator input at Mach number 0.90. 
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Time, t, sec 

Figure 5.- Comparison of flight data with time histories 
computed by using the aerodynamic parameters of 
table IV for elevator input at Mach number 0.98. 
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(a) Pitching-moment derivatives. 

Figure 6.- Comparison of extracted derivatives with Mach number for the F-8 supercritical- 
wing aircraft with preliminary results from the Langley 8-foot transonic pressure tunnel. 



O 8 Transonic Tunnel 
O Present study 



Mach no. Mach no. 

(b) Normal-force derivatives. 

Figure 6.- Concluded. 


